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Abstract 

The stability of isotropic cosmological solutions for two-field models in the Bianchi 
I metric is considered. We prove that the sufficient conditions for the Lyapunov sta- 
bility in the Friedmann-Robertson- Walker metric provide the stability with respect 
to anisotropic perturbations in the Bianchi I metric and with respect to the cold dark 
matter energy density fluctuations. Sufficient conditions for the Lyapunov stability of 
the isotropic fixed points of the system of the Einstein equations have been found. We 
use the superpotential method to construct stable kink-type solutions and obtain suffi- 
cient conditions on the superpotential for the Lyapunov stability of the corresponding 
exact solutions. We analyze the stability of isotropic kink-type solutions for string field 
theory inspired cosmological models. 

1 Introduction 

To specify different components of the cosmic fluid one typically uses a phenomenological 
relation p = wq between the pressure (Lagrangian density) p and the energy density g cor- 
responding to each component of the fluid. The function w is called the state parameter. 
Contemporary observations |T| give strong support that in the Universe, the uniformly dis- 
tributed cosmic fluid with negative pressure, the so-called dark energy, currently dominates 
with a state parameter value approximately equal to —1: 

w DE = -1 ±0.2. 
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Strong restrictions on the anisotropy were found using observations [21 [3] , and it was also 
shown that the Universe is spatially flat at large distances. It has been shown in (U |5l |6] 
(see also reviews [7J |8] and references therein), that the recent analysis of the observation 
data indicates that the time dependent state parameter gives a better fit than wde = — 1? 
corresponding to the cosmological constant. It, in particular, gives reasons for the interest 
in models with wde < — 1- The field theory with wde < —1 is also interesting as a possible 
solution of the cosmological singularity problem [91 [TUl [TTj . 

The standard way to obtain an evolving state parameter is to add scalar fields into 
a cosmological model. Two-field models with the crossing of the cosmological constant 
barrier Wde = —1 are known as quintom models and include one phantom scalar field and 
one ordinary scalar field. Quintom models are being actively studied at present time [TJ 
[121 [ISl Ell ESI EE]. The cosmological models with wde < —1 violate the null energy 
condition (NEC), this violation is generally related to the phantom fields appearing. The 
standard quantization of these models leads to instability, which is physically unacceptable. 
In [T3], the theory with Wde < — 1 was interpreted as an approximation in the framework 
of the fundamental theory. Because the fundamental theory must be stable and must admit 
quantization, this instability can be considered an artefact of the approximation. 

We note the problems of quantum instability that arise in effective theories and are related 
to high order derivatives [171 EH]. In [191 EQ] the instability problem has been reduced to 
the problem of such choose of the effective theory parameters that the instability turns out 
to be essential only at times that are not described in the framework of the effective theory 
approximation. In the mathematical language it means that the terms with higher-order 
derivatives can be treated as corrections essential only at small energies below the physical 
cut-off. This approach implies the possibility to construct a UV completion of the theory 
and to consider these effective theories physically acceptable with the presumption that an 
effective theory admits immersion into a fundamental theory. 

In the Friedmann-Robertson- Walker (FRW) metric the NEC violating models can have 
classically stable solutions cosmology. In particular, there are classically stable solutions 
for ghost models with minimal coupling to gravity. Moreover, there exists an attractor 
behavior in a class of the phantom cosmological models [2TL l22l [14] (attractor solutions for 
inhomogeneous cosmological models were considered in [23]). 

The stability of isotropic solutions in the Bianchi models [2H [251 EE] (see also [27] ) 
has been considered in inflationary models (see [281 129] and references therein for details of 
anisotropic slow-roll inflation). Assuming that the energy conditions are satisfied, it has been 
proved that all initially expanding Bianchi models except type IX approach the de Sitter 
space-time [3U] (see also [211 E21 E31 El])- The Wald theorem [30j shows that for space- 
time of Bianchi types I VIII with a positive cosmological constant and matter satisfying 
the dominant and strong energy conditions, solutions which exist globally in the future have 
certain asymptotic properties at t — > oo. 

The standard way to analyse the stability of solutions for Friedmann equations in quintom 
models uses the change of variables [121 EH E5] (see also [S]). In the case of exponential po- 
tentials such a transformation is useful, because it transforms the class of nontrivial solutions 
to fixed points of the new system of equations [12] . We show that the stability conditions 
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for an arbitrary potential that were found in [14J can be obtained without introducing any 
new variables. 

Studying the stability of solutions in the FRW metric, we specify a form of fluctuations. 
It is interesting to know whether these (isotropic) solutions are stable under the deformation 
of the FRW metric to an anisotropic one, for example, to the Bianchi I metric. In comparison 
with general fluctuations we can get an explicit form of fluctuations in the Bianchi I metric, 
which can probably clarify some nontrivial issues of theories with the NEC violation. 

In this paper we consider the stability of isotropic solutions in the Bianchi I metric. 
Interpreting the solutions of the Friedmann equation as isotropic solutions in the Bianchi I 
metric, we include anisotropic perturbations in our consideration. The stability analysis is 
essentially simplified by a suitable choice of variables [251 1271 135] . In this paper we show 
that for an arbitrary potential the stability conditions, obtained in [H], are sufficient for 
stability not only in the FRW metric, but also in the Bianchi I metric. We also analyse the 
stability with respect to small fluctuations of the initial value of the cold dark matter energy 
density. We consider quintom cosmological models, as well as models with two scalar (or 
two phantom scalar) fields. 

The stability of a continuous solution tending to a fixed point implies the stability of 
this fixed point. Using the Lyapunov theorem [361 EZ] we find conditions under which the 
fixed point and the corresponding kink (or lump) solution are stable. For one field models 
the sufficient conditions for stability of isotropic solutions, which tend to fixed points, have 
been obtained in [38] . 

In [TBI [39] the superpotential method has been used to construct quintom models with 
exact solutions. In this paper we get the stability conditions in terms of the superpotential 
and use the superpotential for construction of two-field models with stable exact solutions. 
We also verify the stability of solutions, obtained in the string field theory (SFT) inspired 
quintom models [T3], EH] ■ 

The paper is organized as follows. In Section 2 we consider two-field model with an 
arbitrary potential and we find the sufficient stability conditions in the FRW and Bianchi 

1 metrics. In Section 3 we remind the superpotential method and obtain conditions on the 
superpotential which are sufficient for the stability of exact solutions. In Section 4 we check 
the stability of kink-type solutions in a few SFT inspired cosmological models. In Section 5 
we summarize the results and make a conclusion. In Appendix we present the stability 
conditions on the superpotential in the case of one-field cosmological models. 

2 Sufficient stability conditions 

2.1 The Einstein equations in the Bianchi I metric 

We consider a two-field cosmological model with the following action 
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where the potential V(<f>, £) is a twice continuously different iable function, which can include 
the cosmological constant A, is the Newtonian gravitational constant (87rGjv = 1/M P , 
where Mp is the Planck mass). Each of the fields <ft and £ is either scalar or phantom scalar 
fields in dependence on signs of the constants G\ and C%- 
Let us consider the Bianchi I metric 

ds 2 = - dt 2 + a\(t)dxl + a\{t)dx\ + a\{t)dx\. (2) 

It is convenient to express Oj in terms of new variables a and fa (we use the notation in |35|): 

ai (t) = a(t)e ftW . (3) 

Imposing the constraint 

/3i + /3 2 + /3 3 = 0, (4) 

one has the following relations 

a(t) = ( ai (t)a 2 (t)a 3 (t)) 1/3 , (5) 

Hi = hi/di = H + fa, and H = a/a = -{H x + H 2 + H 3 ), (6) 

where the dot denotes the time derivative. To obtain ([6]) we have used the following conse- 
quence of (0J: 

/3i+/3 2 + /3 3 = 0. (7) 

Note that fa are not components of a vector and, therefore, are not subjected to the 
Einstein summation rule. In the case of the FRW metric all fa are equal to zero and H is 
the Hubble parameter. Following [251 [35] ( see also [2Z]) we introduce the shear 



a 2 ee 2 + (5 2 2 + ft. (8) 
The Einstein equations have the following form: 

3H 2 - l -a 2 = 8ttG n q, (9) 
2H + 3H 2 + iff 2 = - 8nG N p, (10) 

+=*< *=- w *-kw 

1 dV 

£ = C, (12) 



where 



Q = y 2 + % 2 + 0, P = y 2 + % 2 - V(i>, 0- (13) 
For /3j and a 2 we obtain the following equations 

fa = -3Hfa, (14) 
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|M (15) 

Functions H(t), a 2 (t), <f>(t), £(t), and p(t) can be obtained from equations (I9l)- ({T2"I) and 
( IT5|) . If -ff(t) is known then can be trivially obtained from (TH]) . We show in the next 
subsection that functions H(t), and u 2 (t) are very suitable to analyse the stability of 
isotropic solutions in the Bianchi I metric. 

2.2 Sufficient conditions for the Lyapunov stability of a fixed point 

Summing equations and ffHJj) we obtain the following system 

H = -3H 2 + torG N V(<j>,Z), 
1 dV 

^=-^-clW ( 16 ) 
i = c, 

Let the fields <fi and £ tend to finite limits as t — > +oo. System ffTBT) has a fixed point 
y/ = ( H f> 0/> V>/> £/> Cf) if and only if 

Hj = 8 -^V(^ f ), (17) 
^ = 0, (18) 
Cf = 0, (19) 

n = ^/.o)=o, (20) 

*? = ^(^O) = 0, (21) 

All fixed points yf correspond to i/jf = and = 0. We denote the fixed point yf = 
(Hf, 4>f ,i[)f, 0,0) as yf — (Hf, (j)f,ijjf). To analyse the stability of yf we study the stability of 
this fixed point for the corresponding linearized system of equations and use the Lyapunov 
theorem [36"t 137] . In the neighborhood of yf we have 

H(t) = Hf + eh^ + Oie 2 ), (22) 

cj>{t) = ^ f + eMt) + 0(e 2 ), (23) 

m = eMt) + 0(e 2 ), (24) 

^(t) = a + ^i(t) + 0(e 2 ), (25) 

C(t) = eUt) + 0(e 2 ), (26) 

where £ is a small parameter. 
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To first order in e we obtain the following system of equations 

kit) = -6H f ln(t), (27) 

Mt) = Mt), (28) 

Mt) = -3^(^-1(7;^) + ^)), (29) 

ii(t) = £(*), (30) 

Ci(t) = -3fl>Ci(«)-^(^>i(*) + ^iW), (31) 

where 



<9 2 V „ d 2 V „ d 2 V 

V H = 902^' ^ = (^/> 0)' *W = 

Equation ( 127)) has the following solution 



= b Q e~ 6H f\ (32) 

where &o is a constant. For asymptotic stability of the fixed point i/f the function h(t) should 
tend to zero at t — > oo, therefore, the asymptotic stability requires that the condition Hf > 
be satisfied. 

The system of four first order equations (!2"8]) - (!3~T|) can be written as the following system 
of two second order equations 

Mt) + 3H f fa(t) + {V^Mt) + V££i(*)) = 0, (33) 

fi(f) + SH f Ut) + 7^ {V^Mt) + V£Ut)) = 0. (34) 

In the case when VjL = the system of equations (13"3"]) - (13"1|) becomes a system of two 
independent second order equations. The general solution of this system is as follows: 



^(t) = fa + D 2 tj e-l^ if mj = 



ei(t) = £> 3 e V 2 2V ' C V +J D 4 e V 2 aV ' ° 2 7 if 9#J^ 4^§, 



• = (A» + At) e-f^* if 9#J = 4^ 
where Di, D 2 , D3 and -D4 are arbitrary constants. 
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For the asymptotical stability of the considered fixed point the functions 4>\(t) and 
must converge to at t — > oo. We obtain that at V£ t = the sufficient conditions for the 
asymptotical stability are 

H f >0, 5>0, 3±> . (35) 

Let us consider the case V£ ^ 0. 
1. For 

ru2 ,C l C 2 ( 2 ^ V ii\( n rj2 AV U\ i \rll 2 / C]C 2 (V'(c V^' 



the general solution of (l33l)- (|34|) is 



J D 3e -(|^/+lv /A i+ 2 ^)t + j D 4e -(|^/-|v /A i+ 2 v / ^) 



(36) 



(37) 



m = 2C^ X 

Ci^ - c 2 v% - y^j - c 2 ^) 2 + 

+ 2^ X 

here and further Di, D 2 , D 3 , and -D4 are arbitrary constants and 

Vic VL fVl'c VL\ 2 VL 2 

Ai = 9Hf - 2—f- - 2-^, A 2 =U-^ +4-f. (38) 

For asymptotical stability of the considered fixed point these functions must converge 
to 0. As we can see both functions are just linear combinations of exponents to some 
degrees. To satisfy this condition all these degrees must be negative. It is easy to 
obtain that the sufficient conditions for asymptotic stability are 

V'Jc V'L V'c'cV'L V'lc 2 

H f >(} , yr + -yr > 0, - > — — . (39) 

2. In the case 



r „2_CxC 2 ( 2 /^nC72 4 ^'<A i rW/2 , CiC 2 {Vg 



2 
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the inequality A x ^ is valid and the general solution of (13"3"]) - (}3"1]) is 



fait) = D.e-'^f 1 + D 2 e~'i H fH + D 3 e^ H f + ^^ + 



V2 



(40) 



9H 



/ c7 



H f t 



(41) 



It is easy to show that in this case the sufficient conditions for asymptotic stability of 
the considered fixed point coincide with ( )35l) . 



3. In the case = - - ^) 2 , in other words A 2 = 0, and 

Vj ^ ($H 2 f - ^) (Wj - the inequality A x ^ holds and therefore 



Ci 



V" 



+ 



V— C1C2 

+ I D 2 + [ 1 + ^ ^ ^ e -(H f+ VA l)t/ 2 \ 

The sufficient conditions for the asymptotic stability of the considered fixed point are 



H,>0, | + ^>0. 



(42) 



4. In the case 



„ 2 CiC 2 / 9 4V^ 
V£ 2 = 9H 2 - ^ 



16 



f C 2 



9Hj - 



f Ci 











lc 2 





it is easy to check that the equality Ai = is valid automatically and, therefore, the 
solution of (I33l-(|34D is 



Htt 



V" V" 
C 2 C\ 



D x + D 2 t + D 3 t z + D 4 t 6 ) e~? a ' 



f 1 



+ D 4 



V" V" 
62 Ci 



(43) 
(44) 



For the asymptotic stability of the considered fixed point it is sufficient that the in- 
equality Hf > hold. 
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So, we obtained that in the case ^ the general sufficient conditions of asymptotical 
stability of the fixed point yf = (Hf, <fif, tpf) of the system of equations (fT6l) are as follows 

V& V'L VllVL VZ 2 

H f >0, 7 ^ + 7 ?>0, ^£>^L_. (45) 

It is easy to see that the conditions f )35|) obtained in the case = are equivalent to 
the conditions (j4"5j) . if the equality = is substituted in them. 

So, in the general case including all the above particular cases with definite relations 
between the parameters, the sufficient conditions for asymptotic stability of the fixed point 
Df = (Hf,(pf,iljf) of the system of equations f[T6"j) are fT4o"j) . 

Let us add the cold dark matter (CDM) to our model. One-field models with the CDM 
in the Bianchi I metric have been considered in [38] . The generalization for two-field models 
is straightforward, so we point only the most important steps. If one adds into consideration 
the CDM energy density p m , then system (fT6l) should be modified as follows: 



H = 




- 3H 2 + 8nG N ( 


= 












i = 


c, 




c = 










Pm ' 




- 3Hp m . 



(46) 



Let us consider the possible fixed points of system (1461) . From the last equation of this 
system, it follows that at the fixed point we have either Hf = or p m f = 0. Substituting 
(p2I)-d61) and 

p m (t) = Pmf + ep m (t) + 0(e 2 ), (47) 
into the system (|46p . we obtain the following system to first order in e: 



h(t) = -6H f h(t)+8irG N p m (t), (48) 

p m {t) = -3H f p m (t)-3p mf h(t), (49) 

Mt) = Mt), (50) 

Mt) = -BHMt)-±(vhMt) + v!£i(t)), (51) 

ii{t) = Ci{t), (52) 

Ci(*) = -3fl/Ci(*)-^-Wi^iW+^i(*)), (53) 



It is easy to see that the last four equations of this system coincide with the equations 
fl28"l) - (l3~Tl) . Therefore, the case Hf = can not be analysed with the Lyapunov theorem. Let 
us prove that conditions (|45p are sufficient for the stability of fixed points for models with 
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the CDM. First of all, it follows from Hf ^ that p m f = 0. Solving equations (l4"5j) - (l5T)j) . we 
obtain 

p m (t) = b x e^\ h(t) = boe-^ + ^-e-^\ (54) 

where bo and b\ are arbitrary constants. 

We come to the conclusion that if conditions (|45[) hold, then a solution of system of 
equations (|46p is stable. In other words if conditions (|45p are satisfied, then the solution, 
which is stable in the model without the CDM, is stable with respect to the CDM energy 
density fluctuations as well. 

Let us compare the obtained results with the known results. In [14] the stability of 
fixed points in quintom models (C\ = —1, Ci = 1) with an arbitrary potential has been 
considered in the FRW metric. The authors have presented the Einstein equations in the 
form of a dynamical system having introduced the Hubble-normalized variables 



X * = VE5' H = 7M' y = 7%a< (55) 

and considered a system of ordinary differential equations with differentiation with respect 
to the new variable r = log a 3 instead of t. In our paper we have demonstrated that such 
change of variables is not necessary, because the stability conditions can be easily found in 
the initial variables. 

Proceeding from some physical considerations, the authors of [T4] assumed that Hf > 0, 
we proved that the condition Hf < is sufficient for instability of the fixed point. Note that 
the case Hf = 0, or equivalently V(</>/,£/) = 0, has not been analysed both in [TJ] and in 
our paper. In the case Hf > we and the authors of [14] obtained the same conditions on 
the potential, which guarantee the stability of the fixed point. In our paper we have proved 
that the obtained conditions are sufficient for the stability not only in the FRW metric, but 
also in the Bianchi I metric. 



3 Construction of stable solutions via the superpoten- 
tial method 

3.1 The superpotential for two-field models 

Let us consider the superpotential method [H] (see also [42JH3]) for two-field models. We 
consider the FRW metric and assume that the Hubble parameter H(t) is a function (super- 
potential) of <f>(t) and £(t): 

H(t) = W(<P(t),m), (56) 

and that functions <p(t) and are solutions of the following system of two ordinary differ- 
ential equations 

= F(0,O, £ = G(0,£). (57) 
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Therefore, 

6-— F + — G (-—F+—G 

d(p dt; ' d(f) dt; 

and we can rewrite the Friedmann equations as follows: 

SW 2 = St,G n ( ^F 2 + ^G 2 + V 



™ F+ d -K G= - 4 „G !l( C l F° + C 2 G 2 ) 



From (1591) we obtain 



therefore, 



Also, we have 



QW ^T = 87tG n c i F ^i + + — 

ocp \ ocp ocp ocp 



'CzdG _?F\ = ( 1 dW 
r 1 d d(j) d£J \^G N d dcp 1 



F {c 2 d£ d<p) m ^A<kG n C 2 di + G 

If the functions F and G satisfy the following condition: 

dF _ C 2 dG 

then 

dW dW 
— - = - AttGnC, F, — = - AnG N C 2 G. 
ocp a£ 

Note that using condition (1661) and formulas (1671) . one can verify that 

d 2 W _ d 2 W 

d(f)dt; d^dcp 1 ' 

Therefore, to obtain particular solutions of system (l9l)-(fT2l) it suffices to 
relations 

dW dW 

= - 4irG N C 1( j), — = - 4jcG n C 2 ^ 

ocp o£ 

3 _„„ 1 / 1 fdW\ 2 1 fdW^ 2 ' 



V = t-^W 2 - -—^-r I I ^- I + — 



8nG N 32n 2 G 2 N [C^dcf) J C 2 \ d£ 

be satisfied. 
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3.2 Stability conditions in the superpotential method 

Let us obtain conditions on the superpotential W that are equivalent to conditions f )45|) for 
the corresponding potential V. At the fixed point yf = (Hf,(j)f,ipf) we get 

dW dW 
W f = W(4> f ,£ f ) = H f , W rj = _(^ /) e / ) = 0, W J=_(0 / ,e / ) = O. (71) 

It is easy to see that from conditions ( 17B it follows that 

t? = 0, VI = (72) 



and 



V'U = iqcJ^cp (l2CiCtfrG N W f Wh ~ C 2 W^ 2 - C.W'^ 2 ) , (73) 

*S = 16Ci J 27r2 ^ (^c^wywj - c^ 2 - c 2 w% 2 ) , (74) 

The condition iJj > can be rewritten as Wf > 0. 
If WW = 0, then V£ = and conditions are 

(UnGNdWf - Wfa) Wfa > 0, (l2irG N C 2 W f - W&) > 0. (76) 
In the general case conditions (|4"5]1 are written in the following form 

12C 1 C 2 7rG N (C 2 W^ + OW&Wf > C\ {W^f + 2C 1 C 2 (W^) 2 + C?(W£) 2 , (77) 

(UAC 1 C 2 n 2 G 2 N W] - 12C 1 C 2 nG N (C 2 W^ + C 1 W^)W f + W^W^ - (W'^f) x 

x - (^) 2 ) > 0. 

In the case of one-field models the stability conditions on the superpotential are consid- 
ered in Appendix. 
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4 String field theory inspired cosmological models 

4.1 Quintom models with the sixth degree potential 

The interest in cosmological models related to the open string field theory [H] is caused 
by the possibility to get solutions describing transitions from a perturbed vacuum to the 
true vacuum (so-called rolling solutions [IS])- After all massive fields (or some of the lower 
massive fields) are integrated out by means of equations of motion, the open string tachyon 
acquires a potential with a nontrivial vacuum, corresponding to a minimum of the energy 
. The dark energy model [H] (see also [131 HQ [13]) assumes that our Universe is a slowly 
decaying D3-brane and its dynamics is described by the open string tachyon mode. For 
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the Neveu-Schwarz-Ramond (NSR) open fermionic string with the GSO(— ) sector |47J in a 
reasonable approximation, one gets the Mexican hat potential for the tachyon field (see [H] 
for a review). Rolling of the tachyon from the unstable perturbative extremum towards 
this minimum describes, according to the Sen conjecture [48] , the transition of an unstable 
D-brane to a true vacuum. In fact one gets a nonlocal potential with a string scale as 
a parameter of nonlocality. After a suitable field redefinition the potential becomes local, 
meanwhile, the kinetic term becomes non-local. This nonstandard kinetic term has a so- 
called phantomlike behavior and can be approximated by a phantom kinetic term. It has 
been found that the open string tachyon behavior is effectively modelled by a scalar field 
with a negative kinetic term [52]. The back reaction of the brane is determined by dynamics 
of the closed string tachyon. This dynamic can be effectively described by a local scalar field 
£ with an ordinary kinetic term [53] and possibly a nonpolynomial self-action [51] . An exact 
form of the open-closed tachyon interaction is not known. So, following [13], we consider the 
simplest polynomial interaction. 

In the papers [131(35] quintom models {C\ = — C 2 < 0) with effective potentials V(<p, £) 
have been considered The form of these potentials are assumed to be given from the 
SFT within the level truncation scheme. We postulate that the potential is a polynomial. 
Specifically, we assume that the potential V (</>,£) should be an even sixth degree polynomial 

6 6-fc 

v = EE Cfe ^> v ^t) = (79) 

fc=0 j=o 

therefore, if the sum k + j is odd, then Ckj = 0. 

From the SFT we can also assume asymptotic conditions for solutions [T3], [39] . We assume 
that the phantom field 4>(t) smoothly rolls from the unstable perturbative vacuum (0 = 0) 
to a nonperturbative one, for example <fi = A, and stops there. The field £(t) corresponds to 
the closed string and is expected to go asymptotically to zero in the infinite future. Namely, 
we seek such a function <p(t) that 0(0) = and it has a non-zero asymptotic at t — > +oo: 
0(+oo) = A. The function should have zero asymptotic at t — > +oo. In other words, 
we analyse the stability of solutions, tending to a fixed point with <pf = A and £/ = 0. 

4.2 Construction of stable solutions 

Let us assume that there exists a polynomial superpotential W(<p, £), which determines 
potential (179p with formula ( I70j) . To construct an even sixth degree polynomial potential 
V((j), £) we should choose W(<f>,£) as an odd third degree polynomial. Obviously, the suitable 
form of the superpotential is as follows: 

W 3 {<f>, = A>kG n (a 1)O + a 3 , o 3 + a ,if + a ^ 3 + a 2>1 2 £ + a lj2 0£ 2 ) , (80) 
1 Note that quintom models naturally arise from nonlocal cosmological models with quadratic potential |49l 

EBBS]. 
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where a^j are constants. For the superpotential W 3 system (I6"9~j) is as follows: 

(81) 



I 

— (ai >0 + 3a 3jO 2 + 2a 2i i0£ + ai, 2 £ 2 ) 

(^2 



^ = - 7T ( a o.i + 3a , 3 £ 2 + ct2,i0 2 + 2ai )2 0£) . 

Using asymptotic conditions: 0(+oo) = A, £(+oo) = 0, 0(+oo) = £(+oo) = we obtain 

a 10 = - 3a 3fi A 2 , a ,i = - a 2 ,i^ 2 - ( 82 ) 
So, we obtain the following system of equations: 
• 1 



{3a 3 , {<f) 2 -A 2 ) + 2a 2 ^ + a 1 ^ 2 ) , 

^2 

£ = - Ki(0 2 " A 2 ) + 3a , 3 £ 2 + 2a 1)2 



i3) 



and the corresponding superpotential: 

W 3 (0, = AirG N ( - 3a 3 ,o^ 2 + a 3 ,o0 3 - a 2 ,i^ 2 £ + Oo l3 £ 3 + «2,i0 2 £ + oi, 2 ^ 2 ) • (84) 
At the fixed point <fif = A, £/ = 

W f =- 8nG N a 3fi A 3 . (85) 
So, the condition W/ > is equivalent to 03,0^ < 0. It is easy to calculate, that 

^ = 247^03,0^4, W& = 8nG N ai, 2 A, W^ = 87rG N a 2>1 A. (86) 
Using (J73]) - (J75|) . we obtain: 

AA 2 

VU = ~r~ ( 9 ( X - 4nG N C 2 A 2 )al - oJJ , (87) 
^2 

AA 2 

V tt = ~FT («2 1 - <A 2 - l2nG N C 2 A 2 a 3fi a 1>2 ) , (88) 
^2 

= ~~7T~ 2 ~' (3(1 - AitG N C 2 A 2 )a 3 . - a lj2 ) . (89) 

<^2 

If a 2i i = 0, then = and the sufficient conditions for the stability are 

a 3)0 A < 0, 4ttG n C 2 A 2 > 1, a 1)2 (a h2 + 127rG N C 2 a 3fi A 2 ) < 0. (90) 
If a 2,i 7^ 0, then conditions ( 145]) are equivalent to 

a 3i0 A < 0, 

2a 21 - a 2 2 - 12irG N C 2 A 2 a 3fl ai, 2 + 9(A-kG n C 2 A 2 - l)a 2 > 0, 

[3a 30 a 12 - a|J [3(4vrG' A rA 2 C 2 - l)a 30 a 12 - 3Q7rG N A 2 C 2 (4:irG N A 2 C 2 - l)a 2 + a 2 21 ] < 0. 

(91) 

14 



4.3 Examples of stable solutions 

The case of superpotential W(<p, £) with 02,1 = and a ,3 = has been considered in 
In this case the system ( )83l) has the following form 

° 2 * (92) 



x (2cti,2 — 3a 30 ) i 2 , 2a 12 , . 2 tf2 \ / n o\ 



Equation f l9"3"]) can be integrated in quadratures: 

v /e B - 2 (35 + 2)C 2 



rfg = ±(t-to), (94) 

(Y2BA 2 + 8A 2 - 4£ 2 )£ 3B a 2 2 + (35 + 2)C 2 D l 

where Di and to are arbitrary constants and B = (23,0/^1, 2- At = —1/3, a 2; i = 0, and 
a o,3 = 0, the general solution of system fl83|) can be obtained in explicit form [39] : 



A (C 2 e 4a ^ At / c * - Ua\ 2 C 2 A 2 D 2 - 4a 2 A 2 D 2 ) 
M) = — ^ — , (95) 

(C 2e 2a li2 At/C 2 _ W A f + 64jD 2 a 4 £2^2 



IQD^la 2 2 A 2 e 2a ^ At ^ 

Ut) = 5 • (96) 

(C 2e 2ai, 2 At/c 2 _ 2D 2 a lt2 A) 2 + MD 2 al 2 C 2 A 2 

Let us analyse the stability of the exact solution. One can see that <p s (t) and £ s (t) are 
continuous functions, which tend to a fixed point at t — > oo. Therefore, the obtained exact 
solution is attractive if and only if the fixed point is asymptotically stable. At g^o = — Oi, 2 /3 
we obtain that three stability conditions ( 190]) transform into two independent conditions 

a ^ A >0 ' c ^ Scfe?" (97) 

From these conditions it follows that a\ >2 A/C 2 > 0. 

Let us check the stability of solutions, obtained in [39] . In [39] the author has considered 
the quintom model with the following energy density: 



8nGhr ( 1 -o 1 •■> , , \ 



V 

where 

2 = 9o 

p 8nG N 

g Q is the open string coupling constant, M s is the string mass, 

(100) 



^ = 8CT' (99) 
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where w is a nonzero constant. 

Therefore, we obtain C 2 = M| / g a , the potential 




(101) 



corresponds to the superpotential 



W s = 47rG N uj(j) A 




2 



1 



(102) 



This choice corresponds to 



CO 



(103) 



ai,2 



2 A' 



so the stability conditions are 



<jj > 0, 



4ttG n A 2 M* > g, 



(104) 



So, we come to the conclusion that the exact solutions, obtained in [3H] are stable for 
sufficiently large A. 

5 Conclusion 

We have analysed the stability of isotropic solutions for two-field models in the Bianchi I 
metric. Using the Lyapunov theorem we have found sufficient conditions of stability of kink- 
type and lump-type isotropic solutions for two-field models in the Bianchi I metric. The 
obtained results allow us to prove that the exact solutions, found in string inspired phantom 
models [131 EH], are stable. 

Our study of the stability of isotropic solutions for quintom models in the Bianchi I metric 
shows that the NEC is not a necessary condition for classical stability of isotropic solutions. 
In this paper we have shown that the models [131 EH] have stable isotropic solutions and 
that large anisotropy does not appear in these models. It means that considered models 
are acceptable, because they do not violate limits on anisotropic models, obtained from the 
observations [21 E]- 

We also have presented the algorithm for construction of kink- type and lump-type isotropic 
exact stable solutions via the superpotential method. In particular we have formulated the 
stability conditions in terms of superpotential. 

This work is supported in part by state contract of Russian Federal Agency for Science 
and Innovations 02.740.11.5057 and by RFBR grant 08-01-00798. LA. is supported in part 
by grant of the Program for Supporting Leading Scientific Schools NSh-795. 2008.1. S.V. 
is supported in part by grant of the Program for Supporting Leading Scientific Schools 
NSh-1456.2008.2. 
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Appendix. Stability conditions on superpotential in a 
one-field cosmological model 

The main goal of our paper is to consider stable solutions in two-field models. At the same 
time it is convenient to remind the superpotential method for a cosmological model with one 
scalar field <6, which is described with the action 

■S'= / ~ Itt^- ^<r^oi),o---V{o)\. (105) 



IQttGn 2 

where the potential V(<f>) is a twice continuously differentiable function and C is a nonzero 
real constant. 

It has been shown in [3B] that to find sufficient conditions for the stability of the isotropic 
fixed point in the Bianchi I metric one can consider the spatially flat Friedmann-Robertson- 
Worker Universe with 

ds 2 = — dt 2 + d 2 {t) (dx\ + dx\ + dx 2 ^) , 

where a(t) is the scale factor. In the FRW metric the field <fi depends only on time. 
The Friedmann equations can be written in the following form: 

H= - 4nG N C (0) 2 , 3H 2 = 8nG N (f (^) * + V$)j . (106) 

In [38] it has been proven that the fixed point y/ = {Hf, <j>f) is asymptotically stable and, 
therefore, the exact solution (<f>(t),H(t)), which tends to this fixed point, is attractive if: 

— ^ > and H f > 0, (107) 
G 

in this subsection a prime denotes a derivative with respect to </>. 

System of equations H106[) with a polynomial potential V(<f>) is not integrable. At the 
same time it is possible to construct the potential V(4>) and to find H(t) if <p{t) is given 
explicitly. 

Following [H], we assume, that H(t) is a function of <fi(t), called superpotential (for de- 
tails of the Hamilton-Jacobi formulation of the Friedmann equations and the superpotential 

method see also 021 H31 031 ES]), that is H(t) = W{4>(t)). Using equality H = W'ip, where 
W' = we obtain from system ( II 06ft : 

=W', (108) 



4ttG n C 



V = -4^W 2 ~(W') 2 . (109) 

SnG N 32n 2 G 2 N C V J K J 
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The superpotential method is to choose W{4>) in such form that both 0(t) and V(ij>) 
have required properties. Equation (11081) is always solvable at least in quadratures. Formula 
(11091) allows one to find the potential V, provided the superpotential W is given. 

Let 4>(t) tend to a finite limit 0/ at t — > +oo. We assumed that V(4>) is a twice continu- 
ously different iable function, therefore, V(4>f) is finite and Hf = W(4>f) is finite as well. So, 
system (11061) has the fixed point ijf = (Hf, (pf). It is easy to see that 

V'($ f ) = 0, H} = grG N (y($,j). (110) 
From f l 1 9 1) we get the condition 

^ = Z'fn^n (^G N CW(4> f ) - W"(4>f)) = 0- (111) 

If W'(<j)f) 7^ 0, then from f ll08j) it follows that 0/ is not a fixed point, so, we analyse only 
the case W'((f>f) = and obtain that 

= T^n (^G N CW(4>f) ~ ■ (H2) 

16it 2 G z n C v / 

Thus, we come to the conclusion that to construct a stable kink-type solution one should 
find such W(<p) that <p(t) tends to a fixed point 0/ and the following conditions are satisfied 

W"(4> f ) (l2irG N CW(} f ) - W"(i f )\ > and H f = W(} f ) > 0. (113) 

Note that the obtained conditions are sufficient for the stability of the obtained isotropic 
solution in the Bianchi I metric as well as with respect to small fluctuations of the CDM 
energy density [38J. 
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